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CKNKKAI. INTRODUCTION 


TIuMtrct i I'fil fliiiil luci'liunicH rt-Ht’ari'h I al I y Hi*i-kK to fimi 

i‘oni|i lot o qii.iiit i t at i VO and qualitativo iioHcr i |>t iona ol llio volorlty, proa- 
Hurt', timporuuro and iny otiior rolovant prciport y-f lolda of fluid flow 
ayatonui of divorao boundary and initial condltionat and oapocially t«> 
undorat.ind bow tlioao proportyf ioltla dotorroino, nr ari> tbomaol voa inlln- 
onood by, aucb spoclul and oftt>n vury important fluid dynamical phonomona 
aa flow iiiHtabl 1 1 ty • tranaitlon and turbiilonro. 

Karly cfforta in tboorolioul fluid roochanica focuaod »*n tlio ao-oalJi-d 
"potontial flow" of Idoal IncompreaaibU* fiulda. Quito a larno rlaaa ol 
flows could bo doscrlbod aa potontial flowa and tbo analytical niotbuda 
omployod to doacrlbo aucli flows wore almost porfort. Howovor, vlscoaity 
waa sot>n rocounlzed aa a roal fluid dynamical property critical to an 
undorstiuding, .ind to a coraploto doacrlption, of any important flow pbono- 
monon. Tlila knowledge led to the creation of the mathematical model ot 
viscous fluids governed by the basic Navlcr-Stt»kes (N-S) equations. 
Virtually every effort in theoretical fluid mechanics since tlie ftirinula- 
tlon of the N-S equaliona has become one of finding solutions to complete 
tjr simplified versions of the N-S equations for prescribed boundary 
and initial conditions. 

But there in now a rapidly growing belief among fluid nkchanics 
researchers th.il either the N-S equations may not be a completely correct 
malhemutical model o». the general real fluid flow system, or they in.iv 
not have been adequately laideratood by thi-lr users. It appi-ars that 
some corrections need to be ma<le regarillng especially the expressions 
or conception of the flui«l stress-strain relation employed in the N-S 
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iuiih , lliv lumiui.iry I'ouditioiiH ciup I oyt'd for llif v.triabloH, .iiul (tie 
K«-iu‘t>il ii.itiirL' of tiu- varl.thloH t hoDiKc I vom . Kiirthenmiri* , not wl t loU iml- 
the appari-nt incompli-teiu-tiH of clit* doHi rlpt ivo powor, or In (ho 
• indci Hl aiuling «>f tho N-S oquatioiiM <jh u malhoinal ical imxlol ot the ^eiieriil 
fluid I low ayatem, tlie lualheroat leal prohletn of hoIvIiik thoae e<|ti.it iona , 
for any other than Che airoplest caHea, remains a rather formidahie ami 
often yet impoHulhle Uak. From Che few Hitliitionu of the N-S equations 
available, very little intorniation on very important flow phenomena biieh 
aa flow Inatahility, tr<inaltiun and turbulence are yet poaaible. Theat' 
and aimilar rolsgivinKa lead one to the inevitable opinion that it would 
he unlikely to formulate a aatiafactury unified model of general fluid 
flow capable of describing such fluid flinr>f phenomena as turbulence, within 
the framework of deterministic classical mechanics and t liermedynami cs 
such as the ordinary understanding of the N-S equations would represent. 

In otlier words, if the N-S equations are to be employed generally for 
the description of fluid flow systems, it seems that all Che variables 
in the N-S equations must be understood Co be statistical, not determin- 
istic, quantities. 

In the face of these realities theoretical fluid mechanics research 
turned to statistical methods, especially, and often exclusively, li>r 
the description of the turbulence phenomenon. Initially these statis- 
tical methods consisted in direct formulation and application of proba- 
bility distribution functions to describe the turbulent velocity, pre.ssur< 
and temperature fluctuations in simple turbulent flows (refs : 5 , 7 . ii ) • 
More recently, powerful methods (refs: 2., K, 12,, !•♦) based upon statis- 
tical mechanics and thermodynamics have been devi-lopi-d. in which n« w 
iiKidcIs of Che fluid flow system (usually as a set ol interacting | irti- 


cU'h) ari* coni'i-ptual and in wliirh nu'chunlcal t li.ir.ictfrlut it'H .ire 

viewed an |>rt>l>.ib i I i t ieH .iiid their vnlneH appe.ir an m.il liem.il icul expeet- 
at iouH. 

Hut ulthoiiKli these statistical meclianical appro.uhes show definite 
premise toward eventually resolvinK tiie iluld flow prublem, there appe.irs 
to be still a leoA way to the realization of that K^ial. Fundamental 
problems, notoriously those of determiniiiK the proper st resH->st rain rela- 
tion in the fluid flow system and the proper "history" effects to lacl- 
litate closure of the set ot descriptive equ.it ions oi the fluid i low 
system, remain virtually untouched, and continue to palgue every effort, 
however sophisticated, at mathematically model lug the Keneral fluid flow 
system. 

From all practical considerations It seems quite valid to conclude 
that theoretical fluid mechanics research has for a while now been stuck 
ill a dead-end alley and does not appear to possess any sensible exit- 
direction. And the time seems overdue to re-evaluate the progress in 
theoretical fluid mech.anics research, as well as to systematically sel 
up the fluid flow problem towards a complete practical solution. 

The urgent need tor a systematic approach in theoretical fluid 
mechanics research cannot be overemphasized. Hitherto, the t rtuul in tliis 
area of human endeavor has been characterized by what may be described 
metaphorically .as the "band-wagon" mentality: someone at some fioint 

introduces some "luw" method and everyone Jumps into his b.ind-wagon 
with modifications and extensions but with little or no umlerst aiul iiig 
of the fundament.il pliilosophv involved; and when the new method is seen 
to le.id nowhere, a lull appears as the crowd waits for another "n«w" 



nelliod to i*RaT)',o. Sucli an appraoch id rno«pIft«* I y unw.irrantfd and liaN 
lu'fti viTv uxpi'nHlvt* in tcroui ot monotarv ruHt and Individual friiNira- 
tii>iiH, ni>t it> iik-nt ion ltd f rul t leddnvSH in lernid of a Itadlt* uinlt' i nl and- 
lii:{ and dolutitm ui tlie tluid i low prol>li*ni--a problt-tu til i>xt rcim.- Inipor- 
tanci.- tt> nuin'd tct linu logical and guiu-ral dcvolopnienl . 

A aydtemalic appruacli to tluid rovchanicn rcdvarcli Kliould lead lo: 

(a) a better perdpective of the dupremal fluid flow problem, 

(b) the Identit Icat ion of the iniimal probleind fiindanH.>ntal to a 
complete dolutiun ut the dupremal problem, and 

(c) a clear recognition of the poddlhi 1 i I led or ImpoHdlbi I it led- 
of complete dolutiond to the fluid flow problem, ad well ad of 
the nature and type of mathematical looIh ideally Huited to 
tackling the m.ithemat leal modeling oi the fluid flow syatem. 

And the redulLd of duch effortd uhould Infuse definite and needed direclitin 

to theoretical fluid niechanicd research. 

In 19b9 this invest igator undertook this needed Hydtematic investi- 
gation of the general fluid flow problem, and by 1971 hud tully conututed to 
hinidelf that any complete practical uolution of the fluid flow problem 
must Involve the application of a valid and complete "general -system" 
theory. It does indeed appear that the fluid flow system represei'.ts a 
real model of the general dynamical system in nature, to the extent that 
anyone who can effectively analyse the fluid flow system can also, with 
only slight modifications to hid technique, tdtectively describe the 
dynamical characteristics of any othe-r natural system. Following extui- 
slve research in search of a "general-system" theory etimplete and valid 
for all natural systems, this invest igatoi established in 1971 tlu> 
foundation for the formulatit'n of such a theory. Flemints ol this 
"general-system" thet y are presented in reference (1). As will ulti- 
mately be seen, gene ta I -syst era theory either n-solves or points to a 





ilk'tinili- direction for the renolution of nil tlie fundnmeninl problemM 
liitherto flu ountered in llu* complete dcHcription of I Iiild I low KyHlems. 

The present work reprehenis a prellminarv application of lln’ niulir- 
lyin>t principles of this invest iKator's "Keiu ral system*' ilu-ory to tlie 
description and analyses oi tlie fluid flow system. An .illfiiipt is nwide 
herein to establish practical models, or elements thereof, oi the general 
fluid flow system from the point of view of the general system theory 
fundamental principles. Results thus obtained are applied to a simpie 
"(‘xperimenlul fluid flow system," as test case, with parlicul.ir emph.isis 
on the understanding of fluid flow instability, transition and turbulence. 

This report, iiowever, is a presentation only of tlie fundamental 
aspects of the stated work. In later and more detailed ei forts b" this 
investigator, each of the major findings reported herein is taken separ- 
ately and considered in depth. 

The fundamental philosophy that will be employed throughout this 
and future work is that any equations that are to be used must be deiived 
by tibjectlve application of the universal and Invariant general system 
theory to the fluid flow system; and <iny assumptions made must be 
explicitly indicated. Nothing is taken to be sacred unless it conforms 
specif Ically to the universal invariance principles of tlie general system 
theory. That, it seems, is the only way by wtiich we can be sure to 
delect aiul correct all those fundamental errors which hav«- been jirop.i- 
gated tlirough tlie history of fluid iiu-clianics researcli and wlilch iiave 
stymied the development of a true underst ami Ing of fluid I low. 

i'lie "expel Imeiital fluid flow system" cliosen as tin- test case in tliis 
work is the very simple flow satisfying tlie following requirements: 
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(. 1 ) iiu'oRiprt'HHlblf and t iim*-lndc’|>endoitt flaC plate ht>iiiulary 
layer flow gtiveriied by llie equationH: 

du/>tx *■ ;*v/Dy • 0 (l.l) 

ii3u/<)x ♦ v3u/3y • 3(t 3ij/dy)/i«y (1.2) 

(b) I 1 m the relevant ityMtem partlele vlHeOMilv divided by tin- 

maMH denalty: that 1h, If the relevant HyMtem particle 1m the 
fluid molecule, then c • v p/p. It iu aNMUined, with k^mmI 
reaHOit, that the umuuI molecular vlHcuHlty Mhonld nut vary 
In the "experimental fluid flow HyMtem"; 


(c) the following boundary condltionu are MatiHlied. 

u • v • 0 at y ■ 0 . 
u ♦ Uj as y » « 


( 1 . 1 ) 


where u, v, x .ind y have the UMual meanlny,. 
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2 . SI A MoDKi. or Mil. m in n ow svsi i^i 

2.1 liU lOilui l Ion 

WItct I ri-l» r to UH tin* St.it Int leal rj;y A|*( io;n li (M*A) Ih t Ii.u 
Hi'al.’ir appro. ich in tiu* aiwilyMia ol Hyntinnn wliliti Ih Im .r<l upon total 
oiicr^y ari ttu.' primary anJ Indeod hoIo re'ovaiit NyHteni v.iri.ihio .ind wliii'li 
in foiindod upon tIu* tlu'ory of the "mMU-ral-Hyat oi«," an fiiniv latoii by 
thin in 'i'Hl ip,ator in ref t- rone n (1), with HtatlHtlcal nKilianicH ;ih IIm 
primary analytical t<K>l. 

(iuneral Hystem tlieury provldoH the "pliyHical foundation ol tblH 
approach and permitn the extension of a Kciu-ra I i zud mcchanliH to the 
deacrintion of the haulc dynamicH of any Hyntem, once the HyKteni'a 
"relev.int" eleinentH (l.e. HyHtem partlcleH) have been identified. Stat- 
intical ok-chanicH, on the other hand, provides tite analytical tool by 
which the InterartioiiH of the Hyntero particleu and the coimentient Ht.ites 
of the HyHtem iiuiy he studied. 

Details of Htatistical mechanical techniquen for the .in.ilyHin of 
the interactions of seta of particlea are nanally communlv av.iilahle in 
iiKidern textbooks on tliat subject; this inve.st lK*iior especial Iv ti-com- 
meiids the presentations given by Khinchln (U) and Halescu (2.). Hut 
details of the general system theory from which tlic physics of SKA 
derives are beyond the scope of the present work. It .should stiff ici , 
noiu't lieless, to note that the basic tenets of SKA are geiu rally acceptetl 
in nuidern physics .>nd statistics — especially in modern gc leral kinetic 
theory. Sadly, however, the potential value of these tenets in gencr.il 
svsleiii analyses is not very ranch appreciated by ni.iny ol their tisi rs. 

Within the friimework of SKA, the .inslvsis of anv systim revet Is 
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t'tiiivfii it III I V t«* dll' |>toliIi*m ol del vrinln I lift tin* fv<i I ni Ion or itivolotlon 
ol lilt- HVsit-m'H lot.il fiiorKV I iolil hasoti u|>oti .ipp I i iMl> i t< iiiittal •iitil 
tuMiiul.irv roii.H t lonM ; niiy oiitor rflfvnni hyHiom v.ii l.ihW m,»/ i>«> torowfifd 
hv U)>|>ro|»r l.il i>, •nut uHually rolativt iy utM I'lHl.it y tMlrnl .il ioim 

iMsfd upon modal 4tn<ilyHlM of tlio wyMtvni'M tot<il t-nt*i»(V llfid. I'lio h.tHli 
KovorniiiK t‘<]natioit in, ol ('ourwf, tlu* ii>it ol roiihi-rv.it itm ot total 

fiiorgy lor a lyHtom |)iirticlt! liitcrartlnK with Itv "rolovant" oiivl ronro«-nl . 
HowfVvr. it Ih ill tlif luoaiiliiK and the Idi-nt 1 1 Irat ion of the Hyntrm pat ti- 
de lliat SKA Ih moHt Mennitive and ntiHt |>rone to errorn. Nevertlielt hh, 
tile euneeption of llie uyuteni particle can he made exact. 

The model of a Hyntem a«; i. i ronn of interacting Hyntem particieM 
arranp.ed in accordance with Home K^aieral i/.ed Htructiire principle and 
t ranHlorroiiiK In accordance with Home Kcnerall/.ed rot'clianicH Ih what we 
refer to, herein, an the SflA, or Htatintical enerKy, model. I ah.ill now 
briefly prcHcnt the SblA model of the Iluid flow HVHtem. 

2.2 The PhyHlcal Mo del 

KHHent l.il ly , the phyaical SKA model of tlie fluid I low HVHtem in .is 
followH: that the fluid flow uyHtem Ih at anv iiislant a HiatiMtical 

field ol energy Ht.iteH Hpuce-tiroe Htratifled in an order either ol n»in«»- 
tonlcally increaHing or of nitinoltMilcal ly decreaHliiK atate total-eneipy 
valueH. We ni.iv further underHtand thia model when we ronHlder in more 
detail what ft in we refer to hh the enerjty Htate. From tlu* h.isic SKA 
Mkidel ot tlie Keiieral Hyatem, any ayatem la a set ol i nt er.ic t i ok pariicli-a, 
with preacrlhed initial and boundary conditlona. An.l we .lei ine the avateni 
particle (that la, the "relev.int" ayatem elenk.*nt) not neceas.irllv is a 
dlacrete entity hut rather as a atatiatic.il auhap.ice of the auti |ect -ayat eiit 
coni|iriainK all those pliyaical i»r ordin.irv unit a ol the sub Ject -avsi cm 


y. 


whit'll (uiHHfHM |ilt.nt It' ll nuMii IdI.iI fiu-rK'y ninl wIidi - imiI li>ii In th’^cr ilu-ii 
hv i'X.ii't Iv till rt.ii'M' lit .11 I Hi ii'a I i.iw. riif ti’rm "pIivhIi «| nr iirdlii.ii " Ih 
of t'OiirHo, in ri'l.illon lii tin* m.ili* ol |»tri fjit Ion .»t i In* niih )i-t t -Hynifm, 
TIiiih, wliit one .u'ccptM um .i Myi»li>ni'u pnrtic lf ilf|ii'mjH on llio wCiiU* 
wliirli till.* ono |iro|>onfM to doHorthi llii> nyHlom. 

In till* fluid flow ttyHlom tin* phyHicul tir uiiiin.irv uviilfni iiiiitH .ire 
UHUdlly till* I luiti moleciili*«i (oi . In the t'UHe ol tl I Nttor I .it eil IluidH, iloin.s 
or lonixi-il par t 1 1 ' Ii*h) . The I luld I low HyMtem partleli* wliit'h we Mhall 
Himply I'.'ill the "fluid particle" in tliuu any HtatiHtlc.il HubHpace ol the 
fluid flow HyHtem contalnlnti laoleculeH which pohhchh identical me.in toi.il 
eiierKV and wh.me mol Ioiih are ileHcrlbed by exactly the Hame Ht.it iht leal l.iw. 
In other wordH, a fluid particle can be either a single I luld imileciile or 
a group of fluid moleciileH, or In dlHuuclated fluidii, a Hingle atom or 
ionised particle or a group of atomH or lonir.ed partlcleH. 

St ruminHkly (2.0) huH argiM*d quite convincingly that Indeed theoreti- 
cal a» well :iH experimental data indicate that the main difference between 
turbulent flows Involves gcoupa of mole.'uleH rather that lndividu.il mole- 
ciileH; that Ih to Hay, fluid particles In 1 .imi n.~.r I Iowa mav coiiMlKt of 
single moleculi'N while in turbulent flows fluid partlcleH would be gri*ups 
of molecules. Cenerally, of course, the numlier ot fluid niiilK (moleciileH, 
or other) within any fluid particle will depend upon the local total 
energy field. flu* mathematical characterization »il tin- fluid particle 
Hize in a fluid flow system is discussed later. 

Let iiH, however, continue with our physical argument -U ion. Wi- kn.-w, 
for example, that in the absence ot any externally imposi'd energy I iel»l 
Brownian-type motion shiiuld prevail among all the molecules ol .in\ fliilil 
system; that is, all the fluid units would have ldentic.il me.in lot-il 


i‘ii(*ri;y and cXiTiiio Ht .it Int leal 1 v tilMtI.ir wot luit and, tliertlun’ I li> lluld 
iiyatim would lio uiu- larK>* fluid partit'lo. Thun, a xtar.nant I luid xvMti-w 
with no vnort’V input would rtnitain onlv ono iluid (larl li li and tlu loi.oi*, 
would poM«>i‘SH only Internal dynamicH— t he nu>t ion ot tiu u"ita within i ho 
iluid part ii’lo. A Nyntoin can poHMouu external dyinimiiH onlv It It < on- 
tainu nxire than one Hyateni particle. II, however, an external enerp.y 
I ield 1 h ii.ipoHed upon a fluid Myntein, then I low may occur if the ImpoHcd 
enertty lleld la either non-uni form ot 1 m non-uni loriuly louHirained In 
the fluid flow HVMtew: thiH will result in local varl.ilion ot the mean 
total enerp.ies of tSa* fluid units which will manliest as a st i at i I Ic.iinon 
ot the fluid r 'Stem on the hasis of mean total eiieiity — th.it Is, .ei a 
shearinK ot ;ne fluid. In such a flow situ.ition, thoNo Iluid units pus- 
xtrsiiiK identical me. in total enerKV itroup to form the flow energy stateK. 
Within eai h such flow energy state, one can, at least com epi u.il ly , dis- 
tinguish .ind group together those fluid units whose luot ion ohevs I lu' same 
statistical laws--these are the fluid particles. 

Quite ohviously, a fluid particle can only he proh.ih 1 1 1 st I c.i I I y 
described, since by definition it is a subspaie of a iluid space cont.iinini 
those ordinary or physical units (nxilerules in this case) ol tin fluid 
which possess statistically similar motion as well .is Identical iihmii total 
energy. KurihcrnHire, we now h.ive a better understanding ot definite dil 
terences .iiwmg Iluid systems. For Instance, we know that a slagnml 
fluid with no energy input constitutes only one iluiil p.irticle and tluit 
possesses only intern.il dynamics; a sheared I low must he di lined .is one 
wliicli contains .it te.isl tw.i energy slates each ot whiih conl.ilm- Iluid 
particles; and a non-sheared flow must conlain only one eiurgv st iic 
with al least two fluid particles. 


In I'onr iur. ton to ItiiN ptctuMitnl Ion of I Ik- pliyHicnl t oncopt Ion of tin* 
SKA nuuli'l of tfu? fluid flow RyHtom, 1 n-UHt ro-finph.n. I zo that what Ih 
boiiiK dosrrilud 1 m h hta» Iwt icnl , not a doterml n i -,t li , Mold ol onoip.y 
Htutfu. Such u MkMf tat leal field in hvbt concept ua It /od aa an inalant- 
aneoiiH cloud of non-uniform prt»perty dciiMity. ihuu, any p*>lnt in a tiuld 
I K>w HyHtem l» alwaya cnrloHed by n fluid p.irticle; and a fluid particle 
r.innot be conceptualized aa a fixed Hubapace. Thetelore, contrary to 
poHutble critiiiam that the SKA model implien diMcretenean of energy 
ututea in cc itradictiun to accepted coiitiniiuni coiicepta In fluid synteina, 
a careful atudy of the foregoliiK discurtalon ahould obviate the I act that 
the SKA model merely erophaaizea that a m>n-unlform energy lield, if 
imposed upon au undissociated fluid system, fur example, will Introduce a 
stratification of tlie fluid molecules on the liasis >}f nx an total energies 
and tliat us a result of this stratification or "shearing" of the lluld, 
groups of, rather i;ian individual, fluid molecules will beconsf the rele- 
vant characteristic units of the fluid fi/stem. 

2. 1 f’haracter izat loll ot tlie Flow Knergy State 

From foregoing discussions of tlie pliysicul SKA model ol tin- lluld 
flow system, a flow energy state Is seen to be merely .i set of inter- 
acting fluid particles all of whlcli possess identical me. in total energies 
rhuH, one of the characteristic descriptive variables lor anv .arbitrary 
energy sl.ite, p must be the state energy liensity, Cj--ilellned a^: the 
total energy per unit volume in that energy state. The value ol e^ 
sliould be uniform througliout the flow energy stale, p liut , of murse. 
only probabilistically so; for e^ is clearly a stalisiie.il vari.ihle. 

Kach fluid particle within the energy state, p would eoutain inole»’ules 
(or other appropriate physle.il units) whosi* total energy would i’* r i . si'mid 
with e and whose mean total energy would also he iilentical. 


A I low oii«MKy Ht.iti* may liirtln*r lu‘ tioKcr Ibod as a voliiim* of "slniHar" 


IliiiiJ pai I i t' loH . Hv ( lu‘ lorm "Himilat " wi must t'orlainly nusin t liat llio 
tliitd |>arti«'li*H poHsi-HH Idfiitit'al lucaii (ot.il fiii'rp,y imt tli.ir tlu‘ iim)I ions 
ut till* I liii>l particle.! oluy exactly the same staliHliial l.iw. As a lliiid 
Viiliime, a I low energy 8late ronlil he Jescriheil In terim* ol ihara< t er I si i i 
length KcaleH--lor example: its IhickneKs, its length ami its wiiilh or. 

in terms ot corresponding spherical or cylindrical cooidin.ttes. In this 
work, we sImII employ the syinhol , A^, to represent the char.ict er i st i i 
sp.itial size of the arbitrary flow energy state, j*. .ig.iln, A^, most In* 
a Ht.ttistical variahle since a flow energy state is only proiiali I I i st l< a I I *, 
di scr ibahle. 

finally. In a field of Interart ing flt>w energy stati's It wonlil In 
necessary to distingnisli the locat ions of the energy states rel.it ive to 
some felenaice franu*. In this work we shall ciioose as our refermue I r.ime 
the energy state in the field which possesses the lowest energy tlensliy 
value; and we shall call this reference energy slate the "relative zero" 
energy state with J“0. KnrthernHire , we shall employ the symbol, fj, to 
describe the location of any other energy state relative to the zero 
energ\ stale. Moi a precisely, and for compntation.il l.nilitv, the loi.-il 
value of fcj will be defined as the coordinates of any point ot Interest 
in the energy state ^ relative to wh.itever coordinate system hv which tin 
t r.insforiiX'd energy tleld has been described and in which the orit’In ot 
coordinates lies in the lower boundary of the relativ*- /eio eneigv stati. 
in other words, I ^ will be employed as a loc.il po.silion vector r.idl.il to 
t be lower boiind.irv of the relative zero energy state ami eent«-red .it 
the location of tlie flow field being Invest ig.a led. The vector ' ^ most 
.’ig.iin be .1 st.it islieal varlalili* lor .my flow energy state ,. 


Wi* have alri'ady obMt*rvi*d that the fluid f k»w prohU-n coiuiTiis, 
eauvnt f a 1 1 y , tlu* qual 1 tal ivi* and quantitative deacriptiun ol relevant 
property flelda (auch aa velocity, etc.) ol fluid flow ayaleina with 
appropriate elucidation of aaauclated dynamical phen<>nM.-na such us flow 
inatahlllty, transition and turbulence. We seek now to re-forniulate this 
fluid flow problem In the parlance of the statistical energy approach (SKA). 

The SKA model represents a fluid flow system as a I leld of inter- 
acting fli>w energy states; and each flow energy state Is described as a 
field of interacting "similar" fluid particles, where, by the term 
"similar" we mean that the fluid particles possess identical mean total 
energy and their nu)tlons obey exactly the same statistical law. Thus, 
more generally a fluid flow system is modeled by SKA us a field of Inter- 
acting Iluid particles some of which are "similar" and some "oon-slmi 1 ar . " 
Quite clearly this latter conceptualization of the fluid flow system 
readily permits tested imthods of statistical mechanics to be objectively 
.applied in the .analyses of fluid flow. We further noted that a 1 low energy 
state may be characterized by the three ,.tat ist leal variables e^ , Aj, 
and tj , which describe respectively; the total energy density, the spa- 
tial size and the position vector of the energy state. 

Within the framework of the statistical energy approach the fluid 
flow problem therefore reduces to one of describing the statistical 
field fcj for .any fluid flow system, and then recovering any desired pro- 
perty field of the fluid flow system from appropriate secondary calcu- 
lations based upon modal analyses of the energy density field Cj . 

In this work we shall refer to the Instantaneous energy density 
field Cj, as the flow system response energy field. Generally the flow 
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Hynttfin rcHpoiiHo L>nt*rKy fU-ld 1 h a spjci*-t inn* li Ihi r (but ion InvoivinK I be 
vuriabK'H Ap .imi or any other vai i<ibleH iiy wiiirli llif flow eiuTKy atate 
liaH been tii.ir.irti>rt/.L‘J. 

Ttie HtatiHtical enurKy approach thub Hplita the Hiipnmal tluid flow 
problem into two ini inuil problenib, nann'ly: 

(a) Computation of the flow ayatem roaponae eiUTKy iield; and 

(t>) Computation of any dealred Mow Bystom properly Meld hy 
aecondary analyaea ot the reaponae energy field. 

And in each case, statistical mechanica la the baaic and ael f-auKKoat iiiK 

muthem.iMcaI tool. 

Finally, we note that the concept of llie fluid particle introduced 
by SKA roust neceaaarily alter out uaual formulatiuna of transport coeffi** 
cienta in fluid flow aystema. For instance, in tlie transport of moment uin 
we may apeak of irolucular viscosity, per se, only in cases in which there 
exists only internal dynamics, that is, when there exists only one very 
large fluid particle; otherwise, we must define* and employ a "particle" »>r 
"eddy" viscosity and not molecular viscosity. It is clear that such a 
particle or eddy viscosity would always be apace-dependent . Kven in the 
so-called laminar flows the appropriate viscosity la not necessarily tiu* 
usual molecular viscosity, although it seems, from practical experimental 
results, that the appropriate laminar flow particle viscosity equals the 
molecular viscosity in magnitude and distribution. 


2.5 Comparat i V£^ jlem^.jrk^ Abj)ut mid JJther StaJ.istic.il Appr oadies 

At this point I must emphasize that the statistical energy .approach 
described in this work is quite different from current and perhaps more 
familiar statistical mechanical approaches to turbulence studies, siu h 
as are described in references (li, 8, 12, 1^, W, .and 15). Tlu*re .ire 
three major differences that are easily discerned. 
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FiiMt, tlii‘ liaMic rMtilclH of tliv fluid flow HyMti-in differ iHtween SKA 
and ei>tivi'iU ioiia I Hlatiatlial met'liaiiica I methodH. M«)k( , if not all, i on- 
temporary Hlat 1st leal met'haiileal analyneM ot flui<l (lew model the Mold 
flow ByNleiu au a net of interactint> fluid moleiaileH; SKA inodelH the I luld 
flow Hybtem a« a net of interaet ing "groups" ol fluid moU«Miles--lhe fliil«l 
particles — some of which are slnillar and s»)me ol which are m*n-s irol I ar . 

I’t* this investigator's knowledge, only a recent work hy Struminskiy (.to) 
and perhaps somt* other suhsequent work hy the same .ind possilily lt»llow«*r- 
authors have recognized tlie need lor tin- concept i*t lluiil particles. 
Furthermore, as a result of this SKA nH)ilol , tlu* comept ol moleculat vjs- 
cosity gives way to that of a space-time dep«Mulent "pat t li le-v I s<os 1 1 y" ; 
and this realization autonuit ical ly resolves the long-lived stress-strain 
relation problem in fluid dynamics. 

Secondly, the techniques used to simplify the mal heinat ica I analyses 
differ between this work and most other statistical mechanical an.ilyses 
ot the fluid fl»>w system. In this work prohabllltv thet)ry and the methods 
of stochastic au.ilyses «re explicitly emphasized in con t r. ad 1st i net ion to 
the general and special kinetic theory approaches and approximations 
employed by contemporary models. 

Thirdly, and most importantly, unlike most, if not all, other statis- 
tical raerh.inlcal .analyses of the fluid 1 low system, SKA is a strictly 
scalar energy method. The experiences of many invest ig.it ors ovt r many 
ages ot men have shown scalar energy methods to be mote general and more 
powerful than those metlusls of analysis based upon vector concepts of 
force, momentum .and accelerat ion. Nonetheless, 1 nvest I g.it ors had hltheilo 
very reluctantly avoided the use of energy methods primarily bi'i'ause, oven 
though the methods were simple, they gave only global r.-snits .uul in the 
absence of general techniques for the deconiposi I Ion of eiu rgy into otlu r 


ri'ii’VJiil proptTl li‘S HiK li .IS Vfloilly, .a 1 1 > li-r.it Ion , Inrn , loiiipi-i uiiro, 

• U'lisity .111(1 so on, Ihi'sc ciu'rp.y iiu'tluxis conlil mil lx- cniployod in .my 
liol.iiltd woik. Howi'Vft , I ia* gciuT.il sysloni llicoiy i inim i .il cd |iy llil.s 
invoKi iK'it nr has now ii'iiiuvad , or, 11 yon will, poinicil lu .1 dclinilf 
direction tor the removal ul,the arorenienl ioned dll lit nit v; it is now 
possible to recover roost desired property-1 ields I rom .1 nlven tol.il eiierj;y- 
tield. ilniH, the scalar energy met IkkI h.is been eni.tnc 1 p.itt'd and h.is 
regained its 1 ormidabi 1 1 1 y . 

The .statistical energy approach has been applied in dittei'enl torni.-., 
and with reroarkuhle success, to problems in .icoiist ics .aid structural* 
analysis ( U>) . In tact , a numiier of large system dynamical an.ilyses ( l‘>) 
.ire shitting towards SKA especially because of the mathematlc.il siro)ilicitv 
of the method. This investigator strongly anticipates that SfJV will 
become a standard system's analytical method on< e more invest i g. it ors assure 
themselves of the validity yet simplicity of the modal analyses ol the 
system response energy field suggested in this work for the recovery ol 
desired property fields of systems. 

Finally, .SKA very clearly stands higher in the hierarchy ol descrip- 
tive methods than hitherto employed methods ol fluid I low .inalysis and .is 
such, offers a more likely avenue, than other contemporary methods, to 
a better understanding and containment of the turbnlemc- problem .is well 
as (It the problems of flow instability and t r. nisi t ion. To lllustr.ite 
this reality we shall test the SKA model bv applying It to onr sim)ile, 
though, quite practical, "experimental Iluid flow system." 
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i. THK T UKBULE NCK PAf^MFTKRj^ f 

Sl-IA, ruturcMice (I), idcMitlfluH ttu' turbulence tlelil of uiiy syutcro .lu 
tlie net pure fluctuation field of that uyHtem, characterized by a mean 
turbulence energy, e'. And in "generul-Hyatem" nnalyneH the mean value 
or niathemat icnl expectation, e', of tlie turbulence energy field, in hIiowu 
to be a computable fraction, <^(D), of the system'a total kinetic energy, 
e. Ubviuuuly this phi-parameter, which we shall generally call the 
turbulence parameter, is extremely important in the SEA <lescriptlon of 
fluid flow turbulence. We shall investigate In this chapter, the general 
nature of this turbulence parameter, (^(6), with respect to the fluid Now 
system. First, we list the relations and constraints defining |(6), 
as derived in reference (1). 

e ' ■ ♦ (6) • e 

4.(0) » 20exp(-O.7J602) ; and (l.l) 

0 < 4 ( 0 ) < 1 

In the above relations, e' is the local mean turbulence energy, e is 
the local total kinetic energy and 6 is a dimensionless quantity whicli 
can readily be shown to be directly proportional to the local flow fluc- 
tuation Reynolds number. 

From relations (3.1) we deduce, as illustrated in figure (2), that 
4(0) has a maximum value equal to unity when 4 ° 0.U2, and that as 
0 *• 0 or as 0 4^4(6) * 0. Thus, if in any Now system the 0-value 

lies close to unity, most, if not all, of that system's total kinetic 
energy will reside in the turbulence mode. For a standing fluid we 
therefore expect that 0 will assume its critical value of 0.82 and should 
be space- time independent. On the other hand. In fluid flow svsrems we 
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L'xpi-i t tli.it: 

0 < 0 < 0.82; or 

« > I 

Since Wf know that In practical I lultl flow HyMtenin the unit fluctuation 
KeynoldH number 1 h uHually greater than unity and aince 0 in directly 
proportional to the fluctu.itlon KeynoldH number, we m.iy speculate that in 
real flow Kyatema the conatraint, <* '* I, ahoul.l hold. 

3* I I'li* Hhi -Equat ion 

Since the total kinetic energy, e, of a fluid flow avatem ia simply 
the sum of the nH.*an kinetic energy e, and the turbulent kinetic enerj^y 
e', the relation hetwi'en e' and i* , (3.1), yields the tollowing relation 
between e and e. 

e ■ (l-.>)e (3.3) 

We may now establish the equation of evolution of ik(0) by subtracting 
the equation of transport of the mc*an kinetic energy, e, from the equa- 
tion of transport of the total kinetic energy, e, using the relation 
(3.3) above. 

For our simple "experlnknit a 1 fluid flow system" it c.in readily be 
shown that the transport equation for total kinetic energy could be 
written as follows: 

u3e/3x + v3e/9y ■ e/3y^- ( i / 2y-3c/3y) 3e/3y (3.4) 

Substituting tor e .= e/d-if) in equation (3.4) and subtracting the mean 
kinetic equation, as ia customarily done in fluid mechanics research, we 
obtain the equation of evolution of i>(0) as tollowa; 

u3.^/3x + v34'/3y ■ t3^.k/jy^ + {(/(l-'k)) (3<’/3y)^ - !c/2y - 3t./3y (3.^) 

- (c/e) 3e/3y) 3«(./3v 

liquation (I.*)) corresponds exactly to the classical "prey-predator" 
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fqii.it lull i'uiiimunly umplovcd in llu* mIimIv ui ucoHyntuiPN iinii tliu Htiihillty of 
natural HVMt>.ti*H. Thin hIiuiiIiI Iu* i-X|iu('ti‘(i , aiiu'o ^ may bi* vifwoJ as tiu- 
norm.il I zi-it prcilator (tir t urbuK'iici.') pupul.it ion louiliiiK uptm tlio pruy 
(or im>an flow), aiui tlivroby critically dcpundent on tbo qu.intity and qual- 
ity ut the mean flow. 

Tbe folIowin>t bouiul.iry coiuiitioiiH mimt bold: 

♦(0) ► I an y *0 


(3.1,) 


t(0) ► 0 as y » » 


J.2 .‘’J. k**>' 

l.ol iiH now invuHt iKatf Honx- of tlic mon* obvioiiu inipl ic.it ioiiM of tin- 
evolution equ.ition lor tin- turbuli-iu-e p.irnux-te r , , and hi-v to wliat extent 

Huoli impl ie.it ionn conform wit li or contradict expet iim nlal obnervationH 
of f’.uid flow uynteniH. 

We sliall rewrite equation (1.5) In a t r.insf ornx-d form uKln^ t lie i.evy- 

l.eeB t ranHform.it ion for a tlow over curved bodien as fully diucuHHed and 

employed In relereiice (18); and w*> shall uHe exactly tlie notations of 

reference (18). Tluis, the (x, y) space transtorms Into tlie iC(x), n(x, \)1 

space, and the (u, v) velocity field tr.instormH into the (F, V) dlnx-n- 

s ion less velocity field; r describes the liody r.idiuM .ind t describes 

o 

the transverse curvature while ) describes tin- flow type (j'*0, 

for planar flow; )*1 , lor tixlsynuiu-l r Ic flow). Uslnp. the aliove transfor- 
ni.>i ion and with primes indicating partial d i f f erent i.il ion with respect 
to n equation (1.5) transforms to the fol lowing: 

+ C,.!.' + Cj, =. 0 (3.7) 

where: Cp «• l/(l-<|)) 

C, - -IT, -F Tp - T, - T„ - T5I 


Co “ K(H/3C) 
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T| • (^.)/l2y(nll r 't*)! 

To - Vv (»» /p)^/(ct 2 *) 

* V i* 

T, - . •/• 

T,, - 2 K’/K 

T«, • (p/p )'/(p/p ) 

K - -2Cv^,F(p^yp)2/(tl2J) 

l.et UH now (U'paratoly liuipect tlii* cocM I icioiitH Cy, ( anil (Jo oriiir- 
rlntt In equation (' 3 . 7 ). 

C2 in (l>t* cool f iclont ot tlie Hource of non - 1 Iiuar 1 1 y in equation 
( 3 . 7 )J it iK a tuiu'tion only of and iu deiinitely non-zero; tliuH^llie 
evolution of ^(0) t raiiHverse to the flow direetlon will alwavH he non- 
I inear. 

Cj dooB not explicitly contain i, hut 1 h coinpoMed ol live terim. which 
Heem to elfectively import the Influences of various I low conditions as 
fol lows : 

Ti : I'he basic term; also importing the effect of external ly-lnipressed 

Htre.ini curvature 

Tj ; imports the impact of the transverse I low (I.e. transpiration at 
tile boundaries, etc.) and the influence of externally-impressed 
stream curvature; also the direct influence of viscosity and com- 
pressihlllty are carried by this term 

T3 : imports tlie direct effects of shear and thermal stratification 

T4 : imports the effect of slu-ar; and 

: imports the effects of compressibility and thermal stratification. 

Thus, in terms of the Impact or influence ol flow type and I low houndarv 

and initial conditions on the turbulence field ol a flow system, it is 

clearly the coefficient C) that is critical. 

C contains essentially the iniluiuce ot the stre.imwise gradient ol | . 
o 
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Wi< K.iy coiu'IihJo from tla* .ilutvt* tli.il flu* (•voliilion of ^ fntriratoly 
Involvi-H tfio impaclH of flow typo and flow boundary and fnftial cundltloiiM, 
aa wtuild l>o oxpoi'ti-d In real flow HitiialfonK. Howovor, ft rora.ilnM to bo 
dvtormlnod if tbv impai-tH of flt>w type and flow condftfona on tbo turbu- 
loncu paraiiii>ter follow tbu dfroctlona uuually obaorvi-d In roal flow alt- 
ualfona. To maki* Ibia determination we imiat attempt an analytical aolntitm 
of equation (3.7) and then inapeet aueb a aolution, or we may numerically 
experiment on equation (3.7) in the computer, varyintt the different 
Impacta imported by the coefficienta C 2 . C|, and Ci). To an extent we 
aliall execute both alternativea in tbia report. 

An analytical aolution of equation (3.7) may be initiated by trana- 
forminK the equation into the rlaaairal Abel'a equation, with functional 
coefficienta, and employing Kamku'a, reference (10), aolution to the 
Abel equation. Tbeae efforts yield the folU>wing integri>-d i f ferent i al 
equat ion: 


(3.8) 


« -^Ci<)( 1 - 41/2 - (f)‘/ 6 -*-) + (l-(|>)sj 

where: i - /K I /(!-♦)) 11 / (C, f. ) Idf, 

The aolution may be completed for If S can be explicitly inte- 

grated. Nonetlieleaa aome deduct ior.a may yet be made concerning the evo- 
lutions of 4 by the inapectien of tlie partial aolution given in equation 


(3.8). 


We clearly aee from equation (3.8) that the «>-protlle at any ^“Statfon 
will be described by a function of n decreasing fion. a maximum at fixed 
boundaries to a minimum in the free stream: tlie slope o! the curve being 
critically determined by the coefficient Cj . 

The ^“evolution of the turbulence parameter is more difllcult to 
explicitly determine from equation (3.8), but the general nature of the 
solutions to the "prey-predator" equation, to which class equation (3.7) 
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'.li.it f will *>volv4’ likv *' I oj, I Mt 1 1 *' (Mirvi’H, I rni« Inl- 
li.ill) low vjIiu'h throiiKli n t r.iiiHi ( ionni r.iii|..« »>l Kr<rwlii»r, v.iIiu-h to a 
final raat'.tf of htxli, but (’oiiHtaiil or only hIowIv i'lrii;t>l»K> valuou. TIiIh 
evolutionary trend rontornui with phyMlcal reality. 

Kurt hertikirv , by InMpet’tinK the kIkiih oI the t»'rm» 1 1 tlirough 1^ ol 
the coeffltlent C| In the partial Holutlon, equ.ilion (».H), it U easy 
to wee that the direetions of the impaitw ol flow type and I low t'ondlflon.H 
also I'onlorm with phywieal re.ility. 

Thiiw, the turbulenre paranuter, appearw to lie a very plauwlble 
concept ion. 

J. 3 Num erical Solution of the IMil-Kq^uat Ion 

A prel linln.’iry nuim.*rlcal wolutlon ol the t r.inMli>ruu-d equation (1.7) 
iw attempted herein aw followw: 

Uwlng the notation In relerence (18), the lollowlnp approximal lonw are ma.l»'. 
nr+l, n ^*^nrfl, ii+l ^^'^iirKl, n ^'^nrfl, n-l 

^ niKl, n " n+ 1 ^^^nrfl, n n-l 

(!,' )^ i* 

iirH , n ^ ni , n m+ 1 . n 

' » J, ^ . 

nH-l, n m, n • nrfl , n 
(2Af,2 • X;,) 

liquation (1.7) can then he rewritten ai. U>lloww: 

^n ^nrH , n+l ^ **n '^m+1, n ^'n ^iiH-1. n-l m, ii 

A 

wluTe: A ■ -(Yi + (! Y|, ) 

n 

- I Co* (X,Xs - X,)/(-’X.,A.:2) - Y.. - C*Yil 

t:_^ - (G*Y6 - Y,) 


( 1 . 9 ) 

( 3 . 10 ) 

( 3 . 11 ) 

( 3 . 12 ) 

( 3 . 13 ) 

( ». U) 


M, II 


2 1 . 


>> - t-'o (X X., - X,X^)♦ /(2XsAr..,. 

n^“-V .. 

Co’ - -■2tv,Vp)^K/(...')l^,_ „ 

Kqu^itiiin ().I4) haH tht* t ridi.iKi>nal matrix form and can readily be Holved 
to yield 

Tile numerical .i|(|>ruximat Iona empluved lifr<*ln are clearly not the 
heiit poMbible; the converKence obtained hax not been aut lafactorv. None- 
tbeltHH, fiKureH (3) and (4) abow a tvpical evolution of ^ in a ahear 
flow. A more detailed numerical analyala of equation (3.7) will be 
performetl in a future work. 
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4 . I riiKBiJi.Kw.K J.i «»w„f i ki.d 

4. 1 TIu* Mi'.in iurbii! flirt- Knc-rgy Klfl<l 

Ah alrt-.itly huliratt-J in tlif iirt-ffii i iik t li.ipicr , ''ni iuTal-HyHttW 
(K-m) litt-orv Kivt-H tlio nuitlu-nbil Iral fxpfrtat Ion, t-' ol Ilit* t urbulfin t* 

• lutrKy lit'lij of any HyMtem aa a I'onputabU* irartlon, of Miat HyNluni'a 
ttifal kiiiflic fiitTKy fii‘N, o. Siiu'f tlu* lot.il kiiiitU' t’lUTKy llotil ol 
a HyatoRi la (Ifflnod In g-a tiu-ory aa tin* aiini ol tiu- t urluilviuo kiiU'tlr 
i-iioiKy, o', and tlio moan flow onorgy c, ttio lolltiwinp, rolationa roaiill: 

o ' ■ • o ■ Y • o (4.1) 

t - ♦/(!-♦) (4.-2) 

Tho abovo rolationa Imply itnit oiioo tlio ti>tal kinoiii- or tlio moan 
onor>;y fiold of a fluid flow ayalom la kn«>wn, tho rorroapoiid inp, mt*an tur- 
biilonoo kinotlc onorKV fiold roav bo oonaldorod known, if tho liold of tho 
tiirbulonco paramo tor, la alao known. 

laually, liow>>vor, wo nuiy know tho total volooity liold or tho na-.in 
voloi'ity of a flow ayatom and wo wi>uH w.ini to oxplioitly oompulo tho 
moan t'irbiilont volooity fiold in aomo pivon ooordlnato iliroctlona. In 

ordor to acrompllali thla task wo would havo to oxploit oquationa (4.1) 

and (4.2) in aeurch of praotio.il and valid f ormu I .it Iona by which wo m.iy 
rooovor volooity fiolda from total kinolio i-noip.y fiolda in .mv doair.iMo 
otiordin.ito dlroolion. hot ua iniw altompt that t.iak for tlio moan lurbiiUiii 
volooity fiold. 

hot ua roproaont tho tot.il, mo. in iiul turbulonoo volooity liolda ol .i 
tliild flow ayatom by tho vootora V. V, ami V', roapoo 1 1 vo I y . In ihroo- 
dimo-naional oartoalan ooordinato ayatom those volooitv fiolda bo» 

V i (u, V, w) 

V (u, V, w) 

V = (o’, v', w') 


( 4 . ») 
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Ill tfriON ot I lifMc vt'loi'iiy fU'liiM, the- ciuTgloN , i- , o, ;iiul i*' nwiy lu- 
.iH IoIIiiwh: 


4* 

|) ( U + V 

4- 

w )/2 
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i p(u + V 

■f 

w )/2 
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o' 

p(u' + V* 


4 w’ 


(4.4) 


And from itiiiH (A.l), (4.Z). mul (4.H) we obtain the relations: 

2 2 2 2 2 2 >2 
(ii* + V* + w' ) ■ ^iu + V + w i = i <C/ ( 1~^) H II + V + w ) (4.^) 

A general expresHlon for tlie turbulence components in llie principal 

directions of the reference frame mtiy be ^iven as follows: 


.'2 


J 

where; 


•• ?l: 1. 1. - X, y, z 


l)“i 

j 

X 

V, u 


IJ i 


(4.f>) 


u , u 


V. »• 


w; u 


• . 


u , u 


w : u 


U, 


y ■ z X ■ y ■ z X 

w; and the vtcTinH riip resign t the fractional redistri- 
bution of the turbulence energy from one principal coordinate direct loa 
to another. Obviously, the following constraint applies: 


(v + v fy ) ■ (i + Y Y )“(y +Y ■*’Y ) 

’xx ’xy xz yx yy yz zx zy z/. 


(4.7) 


As a first approximation In this work we had assumed that: 


Y, ,• 1 for 1 - 1 

J : 1, J - X, y, z (4. Ha) 

Yjj - 0 for i ^ j 

W»* have found this assumption not to be acceptable since we know that 
w'^ may be quite different from zero even In flows which are two- 
dimensional in the mean. Wt could correct for such an assumption by 
complicating the i(>(6)-functlon; but that would be highly undesirable. 

The Y~terms may be formulated in terms of the total and mean velo- 
city fields by invoking an analogy between the v-tensor, apparent I rom 
the above relation, and the well-known strain tensor. This exercise is. 
however, rather delicate; It will be given full attention in a later work. 
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K>r till* piir|>ONi> of tlu- roino tiuior t'** tliiH work we nIi.iII m.iko •! him-oiuI 
opproxim.it f on, natni'ly, that: 

Y . ■ a 1 ; lor all I ■ | 

; I . ). • X, y, * (A.hb) 

Y j I ■ (l-a)/2; I t)r all I J 

TIiIm aHHunipt ion automat iciil ly HatirtllL-K the t'onbtraintH (4.7) .imJ appearn, 
in fact, to he a reaHonable approximation In nianv prartieal fluid flow 
Hyntems. The quantity, a, appearH to he a syntem eoiiHlant. V.'e hImI I 
experinunt with a few values of a, in tlie r inp,e 0.3 _ a - O.h; hut in a 
later and more detailed work we shall attempt to find appropriate ior- 
mulations for a. 

The foregoing results indiiate tliat once the total or the mean 
velocity field of a fluid flow system is known, then only a knowledge 
of the <f~parameter is required to obtain tlie corresponding mean turbulence 
field for the system. 

4.2 Ij^e I’rohahf lil y Law of Turbulence Intensity 

Oeneral-system theory indicates that if we define a turbulence 
Intensity, I ■ e'/e, the local ratio of tlie turliulence kinetic energy to 
tlie total kinetic energy, tlu'ii tlie proliahllity density lunctlon ol I is 
given as follows: 

I (i) - f (g,)ldg,/di| + f„(g.) |dg^ /di I , 0 < I < 1 
‘ ^ * 14.‘») 


■ 0 ; elsewlleri- 

where: gj ■ iM I + (l-i*) ^}/i 

g2 on - (l-i^)'*'}/! 

f 4») - 128x^/lexp(5 x) -1 f ; o < x < “< 

A 

■ 0 ; elsewliere 


and 0 Ik, as previously defined, proportional to tiie local fluctuation 


Koyn«>ldH iiiimlicr. 

Sliu'i* > “ 2t»cxp(-(). 7360^), our oarlier Hpoiulat Ion tliat tht- cunutraint 
0 < I, niuHl hold /ieldH, tliat : 

0 » {-tn(O.b^)}* (4.1U) 

Thus, oncu wo have computed the ^-fleld fur a I luid flow system we 
nuiy not only compute the mean turbulence field but also, from knowledge 
ot 1^(1), we may completely describe tlie statistical I ield of the turbu- 
lence. 

^ ^ Turbulen ce En ert;y 

Finally, general-system tlieory also provides a spectral representa- 
tion of the turbulence energy field. 

11 in any energy state of a flow system the n¥)st probable turbulence 
frequency is w , then the turbulence kinetic energy, E , which resides in 
the turbulence frequency w " z • in that energy state is given by 
general system tlieory as follows: 


(i:^/e) - 256Hz*’/ Qz^ + f>^){exp(5z) - Uj 

(4.1 la) 

or (E^/e) » 256Hz*’/ ^ (1-4)) (z^ + 0^Hexp(5z) - 1^ 

Employing equation (4.10) we have that: 

(4.12a) 

(Kj/e) =» 256z‘’/D 

(4.1 lb) 

(i:i/3) » 256z*’/{l)(l-4.)l 

(4. I?b) 


where: I> \ (z? - tnO. 54>) / ^fnO. l'*)! expfiz) - 1} 


Quite obviously, these spectral forms conform with observed reality; the 
bulk of the turbulence kinetic energy is seen ti> reside in the lower 


frequency f luctuat ional inodes. 



OH 


■i. nil KI.OW KiM.i) 


In llii> coiilc'Xl ol "Krn*fr jl-Kystt-m” tlii*ory tin* iinsin iiK>tlun of a 
tliiiil 1 low uvKtfin (l.o. till* inoan flow ticid) is tlio Insl .intsiu-ous sum ul 
till' |)uri' t r.ins lat ional and the- piirf rotational nH>tion ol that systoni. 
This Is the usually ohsotved moan mutlnn ol tin.* fluid I low Hystoin. 

l.ot us again denote the total motii>n of tin- fluid flow system by 
the Velocity field, V (u, V, w) , in three-dimensional cartesian coord- 
inate system; and let us denote the mean (or observed) motion by V 
(u, V, w); the pure rotational motion by the vel»»clty and the pure 

translational motion by the velocity, V,. (d , <1 , q )• 

Q K y z 

The mean flow energy, e, is then given as tolii>ws: 


e r p(V’2 + vJ)/2 

H P(1 -^)V^/2 


(5.1) 


► ^ 

The pure rotational velocity field nujy be written as: ■ mxr, where 

in this case r - Ir • i + r • J + r • kl is the vector ot the local 

X — y 7 . — 

fluid particle rotation. Since the general motion t)f a system partli'le 

within any energy state may usually be replaced by the motion of the 

center of voluims A, ot that particle plus a pure rotation and a pure 

vibration of the* particle about A such that the* total kinetic em*rgy 

density of the particle corresponJs to the total kinetic energy density 

value in the e.tergy state, it becomes obvious tlial in this case r may be 

* 

considered identical with the fluid particle characteristic tadiu£ , r^ , 

alieady discussed in reference (>); |i^.| essentl.il ly represents the 

characteristic radius of that fluid particle whose surface passes through 

the sp.ice-point of interest. Keference (1) suggests the folli»wing lorimi- 

lat ions for r : 
c 


- m’x*/(Ro |ii*|) 
^ X 


(5.2a) 


r*' ■ ni 'xv/(Ke I V | ) 
y x' ' 


or: r » niyi ( !-4>)u* ) *| au*/Dy | ) 

y ^ 




• nr xz/(Rf^|w I) 


(5.2r) 


* * * 

wlu‘rt> u , V , w art* the component ti>tal velocities non-dimeiiKionni ized 


with U , and Re *• (U x/v); m is a nuiiM.>rical coiiHtaiit and x, y, and /. 
max X max * ’ 


are respectively the elfective Htreaniwise, transverse and lateral coor«ll- 
nates in three-dimensional rectangular cartesian coordinate systi>m. 

The pure rotational velocity may also he written In three-d InM'nslonal 
cartesian coordinate system as follows: 

V„"(uir - u»r)*l + (wr - u»r)*J + (uir - wr)*k (5. J) 

R yz zy — zx xz -* xy yz - 

where: 


i^x ■ (f*w/3y - 3v/Dz)/2 

Wy ■ (3u/3z - Dw/3x)/2 (5.4) 

» (Dv/3x - 3u/3y)/2 


Thus, the observed mean velocity field may be written as follows: 


u"Q +u»r - wr 
X y z z y 


v»q + wr - wr 

V Z X X z 


(5.5) 


w”ci +u»r - wr 
z XV y X 


and : 

(u‘- f V- + w*’ ) * ( 1 - <(>)(ti^ + v^ + w ) (5.(>) 

From equation (5.6) we may generalize the following relati<nis 
between the observed and the total velocity fields in fluid flow systems: 
u ” {(l-<>> ) • u^ - 4 >Y - Av • w' 


iO. 


V • |-|Y “■ > ) 

yx yy 



I-* 


C).7a) 

" “ ^'"'^''zz^ • w^l-' 

wluTo .iiul the |-li;rniH are as previously Ji-Niietl In tlie prere<l I np, chap- 
ter. If, according to our sect)nd approxlmat ion of the prect'dluK cliapter, 
we take ■ a, f»>r all i ■ J and y^j ■ (l-a)/2 for all 1 ^ J, tlun 
equations (b,7n) reduce as follows: 

u ■ Gj{u^ - G 2 (v^ + w' ))^ 

V “ Gj f v^ - G 2 (u^ 4 w^) 1 ^ (5.7h) 

w ■ Gjlw^ - G 2 (u^ 4 v^))' 


where : 


Gj « (1 - (fa)'* and G 2 - 4i(l - a)/{ 2 (l- 4 a)l 
If, therefore, the total velocity field of tlie fluid flow system 
is ktu)wn, tlte observed mean velocity field may readily he computed from 
a knowledt;e of 41 and a. If we want to compute the actual pute transla- 
tional and |)ure rotational velocities we must employ equations (’i. i) and 


(5.4). 


It is easy to see from the fore^oinn that in the To-called boundary 
layer flows, the observed mean velocity field would be mostly rotational 
velocity: tliat is, pure translation is very sm.il 1 in tiie boundary layer 
compared to pure rotation. 

A 1 1 ernat i vi- 1 y , the qu.int it Ics u , v , and w , iiujy be c»»mputed frt>m 
tlie evolutionary equation for e as olitalneil I rom ecpiat ions (Hh) and 
( S O) in U preceding chapter. For our "expi-r imental Iluid flow 
system" we would have the following equations for e and its components: 
ude/i)x + vde/dy “ Cif^e/Dy’ - '^k/2y - c/(l-4) • '>'l/^y - 'i/:>^'^e/>y (5. da) 

u3u^/3x 4 vdu^/3y t3^u‘^/3y^ - ^c/2y - r/(l-4i) • 3(fi3y - 3r. /'»yj 3u' /3y fj. 8b) 

(u" + v ^)/2 


e =» 


(5.8c) 
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THK a)MI’irrATI()N OK THK. .KU'II) FLOW SYSII M 

b.l A of OcHor Ipl i vo j’h^l loHoplij^ 

1lio fluid flow HyHtoro tn.iy hu dcHcriliiil in any om* of sovoral waya 

I 

dopc-mlinK Hcali* of I tindainont 4>1 i ty of |>ropcrly-f ii'ldu that om* is ' 

interoKtod in. Thiii f rovdoro in done ript ion hcIh up what I have callod a 
hlerarcliy of descriptive philosophy, based upon the eletnentar I tv of the | 

cliosen basic property-field. 

In physical reality the Fluid flow svstem is most t '.indaMentul ly J 

an energy system; that is, its sole basic proper tv-f ield is "energy. " ' 

4 

Any description of the fluid fl»iw system as an energy field may therefore 

i 

be referred to as a "single- f ield" or "general system" approach; it would I 

be the highest physical description of the fluid flow system. 

I.ower down in the hierarchy of descriptive philosophy we may chot)se 
to describe a fluid flow system as a mixture of "forces" and "energy"; 
such a description would conform with the classical Hamiltonian mixed- 
field approach. 

But by far the most popular and, therefore, the dominant descrip- 
tion of the fluid flow system in contemporary fluid mechanics res»>irch 
is the Newtonian mixed-field approach in which the systiin is viewed as 
a mixture of 'momenta" and "iorce" llehls. Clearly, I he Newtonian 
approach is lower in the hierarchy of descriptive ptiilosophv than the 
Hamiltonian appro.ich. 

The reason tor a descent in the hierarchy of descriptive philosophy 
is usu.illy the desire to obtain mtiri* detailed inlormallon in an explicit 
form and on as many property-fields ol tlte subject-system as ptissihle. 

For instance, the general system approach would usually yield only the 

i- 

• V. - J 


i2 


total iMiiTKy I tel J of a HiibJecI Hysiem wtii li* the li.iiiii I ton Ian .ippioach 
would provide the forci field jh well an the total kinetic energy llelda; 
the Newtonian approach, on the other hand, would i‘xplicitly yield the 
velocity, ip.!8ii-denM i ty , t emnerat ute , preHNure and other property-f ielda 
ui the Haine Hul>Ject Kyatem. 

Hut what one in explicit 'Metal I" the one loaea in explicit 

"complexity" of the chosen nx-thod. Tin Newtonian mixed-iield approach 
requires more complicated equations, often with severe closure prohleF's, 
for tlie description of the fluid flow system than dot‘s tlie Hami I tor.lan 
mixed-lleld approach; and, as we have seen in re^ ft) ol this leport, 
the generai system appro.ich descrihes the fluid flow system by one very 
simple equation — the equation of conservation of total energy. 

6.2 T he Ge neral-System Approach 

In employing; SKA (or the system approach) )or the computa- 

tion of the fluid flow system, essentially only the lol lowing three 
procedures are required: 

(O transform the ordinary sp.ice of the fluid flow system by 

sotiK' set of transformations into the p,eiu ral system space, 
ior instance as discussed in references (17) and (I''*-) , 

(11) solve the simple fundamental equation ot the general system 
in the transformed or general system space: ami 

(iii) transform the solutions I rum (li), nboVi>, hack into the 
ordinary space of the fluid system. 

What one obtains in (iii), above. Is the response total inergy f U Id in 
the ordinary space of the fluid flow s; si». Iti order to recover any 
other proporty-f leld from this total energy field one must unplov the 
results of the analyses of energy suggested by this 1 iivi-st igator , some 
of whicli have been presented In this rept>rt, aiul all oi which are based 
upon this investigator's general svstem theory. 
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Tlu> r«)ri'KoinK notwltliMtuiulliiK. •> I uinJ.'iini'nl al dlltirulty in tin* lull 
application of the Hyutem approucli Ht I 1 1 reinaiiiH uiiri-Hol ved at 

the tiiiK* oi thlH report. TIiIh coiii'eriiK the el iril.it i»n of a ainiple Net 
of t r JitHformat ionH fur t ran8furnilii)> the ordinary Np.ice of any Hyutem 
into the p.eneral Hyntein ripace. A numher of relevant t ratiHf orm.it lonH 
have been HUKK^'Hted by different invent i^atorH Nurii a.s in referenceH (II) 
and (3.^, hut thene t ranHfurni.it ioiiH Htill nuike nnire dem.indM on computer 
time and Kp.ice than the preaent inventiKator stroiiKly lu 1 ieves in nereHS- 
ary; alter all, the aim of the exercine of the general Hystem appronrh 
1 h to attain maximum simplicity with m.iximum aeeur.iey and not }uHt to 
find .mother alternative approach. At the moment, ttiven a fairly laiKe 
computer space and time (though much stiviiler than what would he required 
by other contemporary methods to provide the same or equivalent aiiK'iint 
of results) this investigator's general system approach can yield very 
detailed results on any fluid flow system; but such computer space and 
time are not easily available, 
f). 3 A Derived^ Newtonia n Approach 

An .ilternative way to test or demonstrate this invest jg.iior's geii~ 
oral system approach is to deduce from it tlie relevant descriptive 
equations in a lower scale of description, su<;h as the popular Newtonl.in 
mlxed-lleld scale, for instance, and then to examine the correspondeiue 
of the deduced description with experimental Iv verifiable knowleilge on 
real fluid flow systems. In doing so wi- have come to the conclusion that 
the system of N-S equations of fluid mecli.inics is in itseli a complete 
and correct description of the dyn.imics of fluid flow In a vector spate, 
provided that: 

(a) all the dependent v.iriables in the N-S equations are under- 
stood to be total variables admitting only of statistic.il 
descript ion; 
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(|i) tlitf viMcoHlty tfrm, c, is iiniifr<iloo(i {■> ho not tho iiKual 

kinom.itio viHrouiCy b.iHoii upon imilooul.ir ttaiiNport i>f iiK>niont um, 
hut ratlior tho gonoralizeU "fluid p.irililo" or "iddy" vIhooh-' 
Ity Patiod u|ion tho tratmport of motm-iit urn hy tlio coiuoptual 
fluid particloM of fluid flow HyHtimH: ntolot'iilar vlstoaity 
would ho only a Hpoclal oaso of iI.Ih oddy vlHcoaity; and 

(0) ’ tho convontlonal undorutandiiig rl prousuro 1 h modlflod tt» 

explicitly dtHtingulHli hotwoon tho internal pioHUuro ot a 
fluid flow syatom and tho internal gravitational ijrcou in 
tho fluid syKtom: thla ia o.spcr tally important in tho gonoral- 
ization of tho Bornoulli equation. 

Kor tho purpoae of practical computatioiiH in thiu report tho N-S 
equal iuiiH will, thoroloro, bo at copied with thi- vlKcoHity term, l, 
given hy the following aimple form: 

l“v(I+T*i^/m) (h.l) 

whore T is proportional to tho integratou value tho <f-paraitk*tor at 
any x-statlon and is any one of tho currently omployod eddy viscosity 
models in turbulent tliild fli>w analyses. In effect, T will ».orve In 
this case as the usual Intornii t toncy factor for a flow in wlilch transi- 
tion is imagined to begin at the loading edge or entrance region. 

Although equat ion(6. I ) is not exact, it slu>uld suffice to prove 
the validity or invalidity of tho concepts suggested by the general 
system tlieory as applied to our "experimental fluid fliiw system" 

To compute tl>e "experimental fluid flow system" we shall execute 
the following procedures: 

(1) solve the following system of coupled equations for the 
velocity and i^-profilcs, for all x-slations In I lie fliiltl 


flow aystera: 

ju/3x + 3v/3y ■ 0 (h.J) 

u<Hi/(4x + v3u/3y ■ 3 (c3u/3y) /f*v fh. J) 

t ■ v(l + T e^/p) (b.4) 

u3<i>/3x + vD(>/3y ■ i 3^t/9y^ + lt/( 1-f) I (a<i/3y)- 

-{c/2y - 3t/3y - (t /e) i»i /'3y 1 3>f /'^y (h.5) 


35 . 


t « li ^ ♦'■'y (b.h) 

•• f w*’ ■ u* (in thin caNL') (6.7n) 

and L in Kiv^‘n liy the uMiial mixinK lein'.tli eddy viNroxity 
model and the f(»liuwinK hotindary eonditiunx are ot>seived: 

♦ ► I 1 

u and V » OJ HH y ♦ 0 

t ♦ 0 and u ► Uj; aa y * “ (6.7b) 

(it) At oath x-stallon compute the mean flow field aa follows: 
u • (^i(u 2 - G 2 V^)'* 

V - Gi(v’ - 

I 

where G j • ( I - • a) (; 2 "^d - u)/{2(I - ♦ • a) I , 

.ind a ia experimentally varied between 0 .) and 0 . 6 . 

(iii) At each x-atatitni compute the turbuK'nt 1 low tield aa 
f ul Iowa : 

u'^ ■ ('13 (n^ + Gi,v' ) 
v'^ - Gj (V^ + C^u*) 
w’^ - G3G^(n^ + v^) 
where G ■ ♦•a, Gi, ■ (1 - a)/2a and a la aa previously defined. 

These three steps complete the computation ol the I low field. 

Other aspects of the flow aystem may be determined in the iiaual manner 
from the mean and the turbulent veli>clty fields. 

Time has not permitted us to include in this report the reaults 
of the actual demonstration of the above computations for our "experi- 
mental fluid flow system." In a later report these results will be 
presented . 


( 6 . 8 ) 

( 6 .'») 


( 6 . 10 ) 
( 6 . 11 ) 
( 6 . 12 ) 


7. CONCLUUINC KHMAKKS 


A prlm^iry conclutiioii from lliiti roHoorrh work in llhit tlitf tliiid I It*w 
uyHtcm Ih Hentilbly, fulJy and practically d»*H«’t Ibfd by tlio 'VA*noral- 
Hyattm" theory enunciated by thia inveat iK^tor . Thia Mlalement renta 
firmly on at leaat the foil«»wlnK thei>retical and practical tibservut itma ; 

Firat, although the explicit derivation ol tlie governinK equationa 
of fluid flow from the fundamental equation of the "genernl-ayateiu*' ia 
not included in this report, tlie aimilarity between the latter and the 

Navier-Stokea (N-S) equationu la rather uhvioua. Indeed, the N-S equ.i- 

» 

tlona can readily be alu>wn to derive from the fundamental c‘quation of 
‘he general ayatem; hence our confidence in the N-S equationu au a valid 
model of fluid dynamica. But, although the N-S equationu may now ep)oy 
our full confidence, we note that without the benefit »>f the inuight into 
their nature generated by "gerieral-uyutem" theory ami tl>e SKA model, 
the N-S equationu may eauily become grouuly mluunderutood and miuuued. 

Secondly, from current theoretical and experimental knowledge ot 
the fluid flow uyutem, tlie SKA model preuented in tliiu work appearu to 
be Inherently conalutent ard complete, vlu-a-viu, the proi^lsion of deep 
insight into the real (ihyulcu of fluids. The nuidel explrinu the n.il 
mechanics of fluid flow In greater d<;tail than do mtist o.'her conti*ropor- 
arv models of the fluid flow uyutem. But above all, by obviating very 
clearly the relationship amcog the diflerent types ami i he ilitferent 
modes of fluid flow, the SKA imidel greatly ulinplilles the analyses ol 
fluid flow uyutemu. In this regard, the naturt* of t iirhu I etice , the 
nature of the "observed" mean t low and the nature of the relat lonuhif* 
bi'tween the turbulence field and the observed (or measuri*d) mean I low 
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fit-hi ilt'Hvt'Vf part I iMil.ir attention. 

lurlnilfiue 1 h clearly identitieU to be tlie puti lltiCtiMtioii.il iNode 
ol the riiiid particleH of the fluid f li>w Myatem. Aiul. contrary to clan- 
Hical allirniationB bucIi an that the turbuJeiice field of the I luhl flow 
MyMtem is not (generally a uniqee function of the mean velocity profile, 
the SKA del demonat raten that Indeed the nn*an tnrhnience field c.m he 
expruHtied .IM .1 computable unique funct ion ol tin* me.in I l«iw I ield. t'eii- 
eral nystera philoHophy clearly expoaea that the errtira In earlier pretlic- 
tion methotlH wliich led to the affirmation that the turbulence I ield 1 h 
not ly ^ uniqii.; function of the mean flow field were not at .til 

due to the reliability of tlu'He met'iods on the mean flow fielil for the 
deHcrlptlon of the turhiilent field; rather, the error« were d<ie to tlie 
fact that iwiHt , If not all currently available theoretical methods In 
fluid mechanicK, are founded either upon no f und.imenta I |>hi losophli al 
cimsiderat iona it all, or up«ni erroneous piii Uisophical wise.icr in>;. Km 
instance, the ex.iminatlon of experimental d.itn should not .lutomat ical I v 
le.id to some intuitive or empiric.il formulation to describe the "ohservetl" 
trend; without consideration of virtually all posslhle variel les ol the 
observed system such intuitlvt* formulations would usual Iv lu‘ .it best 
restrictive .ind would often be conipletelv erroneous, li*r the simple reason 
th.’t dvnaniical system liehavlor is often count iT- 1 nt ul t i ve. To cr«‘at» a 
re.'ilistic relation from physical observations it ’s critical th.it one 
formulate a valid phi losopliical base — a consistent ami complete foiind- 
at Ion or vii‘Wt>oint that inclmles the observed data. Ctmlaiiii'd within or 
tierived from such a phi loso|'hic.il base would be a theory upon which the 
empiric.il or ex.icl relations v.illd In the observed system iiuiy he foimilwi* 
The error-situation beitip, clarified .ihove is .inalop.oos to t ht case 
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ill lilt I Ic iiit' wIkt*' tlii Myniptoms «•! i i .i li.ivt n ll•tll ■..•il ii|mii .iiii| 

ci lit f I y dcHcr Ihcd while the t riu It.iltii< <•( iIm- dtn«a».e .iiui (Ik hh* Ii- 

•init'ii hy wh K'h It )',«'iu*i at eti t hnMi MViiipi nma h.iv« t wi'ipliti ly eliidt'd the 
doctor. Moiit of whut currently nvailahle theoretical ntode1n ut I liild I li>w 
deHt i'ibe are obiiervat lonal— opt leal and inat r iimeiit 1 1-- 1 I Ins I oiih ; hardly 
any model, before SKA, deacribed the real nature of tlie fluid flow Hyatem. 
Thla rather pompoua ataterocuit ia borne out hy the latt (hat, with verv 
little eftorf, virtually every hitherto yiveii deacr ipt ion of fluid I low, 
upon which nuineroua theoretical modela hinKC, can he reconat ruct ed aa 
oliaervat ioiial illualona from the point of view of tlie SKA imid«-l presented 
herein. 

1 would not be belaborinK my criticism of the "aymptumat ic-approa» h" 
commonly practlacd in fluid mechanica reaearch were It not for the Inher- 
ently adverse impacts of that approach. Kspecially, one major line ol 
inquiry, arising from the aforementioned "symptomat Ic-approach ," and om 
which has, more th.in any other, seriously stymied the development of tur- 
bulenc(> theoriea, is the continued combination, in all turbulence iianlela, 
ol the pure fluctuotlon and the |>ure diaaipative energy fields. This 
error and its destructive Imp.icts ar«- not I’asv to see outside the vliw- 
poiiit of the SLA model. It is true that the d i as I p.it i vt> eiurg.y tlel>l 
appears to be derived from the fluctuation energy field, hut it Is also 
true that these two fields are characterist icf I Iv very dlff»-teiit. We nave 
shown in reference (I) that the gross production and gross dissip.ii ion ol 
fluctuation energy in any system ar»* extremely complicated quant it let. to 
theoretically describe. Wli.it we measure as the tiirhulenc'e I iel<l Is oiilv 
the net fluctuation energy field — a much simpler qiiaiiiily to con'i'uli'. 

The dissipative energy field manifests as flu- pressure and tempcratuii 
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ficlilM of any HyKtom and HhuulJ bt* toniilderi’d boparatcly f the lur- 
biiU'nt'o flold: only tlicir inCoructlon or lntord«>|n-ndfn<y nooil be noted. 

Fln.illy, let uu rttiK.-mber that tlie StallaMtal liuTKy Appm^ttb (SKA) 
in a K^'neral approach valid not Juat for the deacrlptlon of the fluid 
flow Hyatem but for the dearription of anv conceivable n.itiiral ayaleiii. It 
la baaed upon an internally cunalatenl and univeraally complete phllo- 
Hophlcal viewpoint t ranacend I tig tlie limitations of sub |ect I v 1 1 y . 
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